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Abstract 



We demonstrate for a generic pseudointegrable billiard that the number of 
periodic orbit families with length less than I increases as nb^l 2 / (a(l)) , where 
60 is a constant and (a(l)} is the average area occupied by these families. 
We also find that (a(l)} increases with / before saturating. Finally, we show 
that periodic orbits provide a good estimate of spectral correlations in the 
corresponding quantum spectrum and thus conclude that diffraction effects 
are not as significant in such studies. 
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Billiards are an interesting and well studied class of hamiltonian systems that display 
a wide variety of dynamical behaviour depending on the shape of the boundary. Of these 
pseudointegrable billiards form a sub-class and these correspond to rational angled polygo- 
nal shaped enclosures with a particle reflecting specularly from the walls They possess 
two constants of motion like their integrable counterparts [|2j but, the invariant surface is 
topologically equivalent to a sphere with multiple holes Q and not a torus. As an example, 
consider the billiard in Fig. 1. For any trajectory, p\ and p 2 are conserved. The invari- 
ant surface consists of four sheets (copies) corresponding to the four possible momenta, 
(±p x ,±p y ) that it can have and the edges of these sheets can be identified such that the 
resulting surface has the topology of a double torus. 

The classical dynamics no longer has the simplicity of an integrable system where a 
transformation to action and angle co-ordinates enables one to solve the global evolution 
equations on a torus. On the other hand, the dynamics is non-chaotic with the only inter- 
esting feature occurring at the singular vertex with internal angle 37r/2. Here, families of 
parallel rays split and traverse different paths, a fact that limits the extent of periodic orbit 
families. This is in contrast to integrable billiards (and to the 7r/2 internal angles in Fig. 1) 
where families of rays do not see the vertex and continue smoothly. 

We shall focus here on the periodic orbits of such systems for they form the central 
object in modern semiclassical theories ||. Not much is however known about the manner 
in which they are organized and the few mathematical results that exist concern the 
asymptotic properties of their proliferation rate. For a sub-class of rational polygons where 
the vertices and edges lie on an integrable lattice (the so called almost-integrable systems 
0]), these asymptotic results are exact. It is known for example that the number of periodic 
orbit families (those which have even number of bounces), N(l), increases quadratically with 
length, I, as I — > oo. For general rational polygons, rigorous results provide bounds on N(l) 
though it is believed that N(l) ~ I 2 even in these cases @||. Numerical verifications of these 
results are few and perhaps none exist for a general rational polygon that does not belong to 
the category of almost-integrable systems. Besides, very little is known about other aspects 
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such as the sum rules obeyed by periodic orbits in contrast to the integrable and chaotic 
limits where these have been well studied 0. 

It is not surprizing then that periodic orbit theories for polygonal billiards have met with 
little success, both in quantizing individual levels ]7|,[8] and in explaining the correlations in 
the quantum spectrum ||. We shall not deal with the question of determining individual 
levels here, but merely point out, that, agreement with the smeared quantum density of 
states has been observed using complex energies |?]|| although alternate convergent schemes 
(such as cycle expansions in chaotic systems are indeed desirable for better resolution of 



individual levels. In this sense, the full scope of periodic orbit quantization in such systems 
is largely unknown, even though several researchers have now looked beyond geometric 
periodic orbit contributions and found evidence of diffractive corrections ||11|| . 

Correlations on the other hand are related to sum rules obeyed by periodic orbits and 
these are robust quantities that do not suffer from acute convergence problems. We thus 
focus on the problem of computing correlations in the quantum spectrum using geometric 
periodic orbits. Further, our results indicate that these estimates are good indicating that 
diffraction effects are not as significant in studies involving the statistical properties of the 
spectrum ||12|| . To appreciate this finding however, it is necessary to study the sum rules 
obeyed by periodic orbit families and show that there are important differences from the 
integrable case contrary to what can be expected from the asymptotic properties of N(l) 
||. To this end, we first verify that periodic orbits in a generic pseudointegrable system 
obey the sum rule [|13|] (J2 P J2^Li a p 8{l — rl p )/{rl p )) = 2nb , where b is a constant, (.) 
denotes the average value, l p is the length and a p the area of a primitive periodic orbit 
family. This establishes the proliferation law, Nil) = nbol 2 / {a(l)) , derived in |13| where 



(a(/)) is the average area occupied by families having lengths rl p < I. Further, we explore 
the behaviour of (a(l)) as a function of length, I, and find that it increases initially before 
saturating to a value much smaller than the maximum possible area spanned by a single 
family. The proliferation law, N(l), is thus quadratic only asymptotically and the number 
of orbits is much larger than that of an equivalent integrable system with the same area. 



For smaller lengths, N(l) is sub- quadratic and this is significant for a two-point correlation 
of the quantum spectrum that we study and for which we demonstrate that periodic orbits 
provide accurate estimates. 

The L-shaped billiard of Fig. (1) that we choose has approximately unit area and has 
no periodic orbit with odd number of bounces at the boundary. It does not belong to the 
class of almost-integrable billiards and is generic in the sense that all sides are irrationally 
related and periodic orbit lengths are non-degenerate. Unlike some degenerate cases @ 
where periodic orbits can be labelled by two integers (analogous to the winding numbers in 
2-dimensional integrable systems) , orbits in generic L-shaped billiards are described by a set 
of four integers though not every point on the 4-dimensional lattice corresponds to a real 
periodic orbit. This makes it difficult to study sum rules from purely classical considerations. 



An alternate approach adopted in [13] uses the semiclassical trace formula which expresses 



the density of quantum energy eigenvalues in terms of periodic orbits JT|: 
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where Jo is a Bessel function, {E n } and d av (E) are the quantum energy eigenvalues and 
their average density respectively, and k = \JE [Q. For convenience, we have chosen the 
mass, m = 1/2 and H — 1. Starting with : 
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9(l) = Y,fW E n,P)= dEf(E;P)Y,$(E-E n ) (2) 
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where f(vE] 0) = Jo{\fEl)z~^ E an d < e < E , it is possible to show using eq. flU) that 
for p -> 0+: 
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where b = E P EZi So dEJ (VEl)J (VErl p ). 

It is argued in ]13[ that bo is a constant and we demonstrate this here by plotting 
S(l) = J2pJ2 r ^t m Fig- 2 where the summation is restricted to all periodic orbits with 
rl p < I. It follows from Eq. © that S(l) ~ 27rb l. 



Fig. 2 shows the behaviour of S(l) for a rectangular (integrable) and an L-shaped (pseu- 
dointegrable) billiard. In the former case, the orbit lengths, l p can be expressed in terms of 



the winding numbers (M,N) on the torus [|T^,|T^] while the areas {a p } are four times the area, 
A, of the billiard except for bouncing ball orbits for which they are twice the area. For the 
pseudointegrable billiard, both the lengths and areas are determined numerically using two 
different methods. We illustrate one of these by first noting that one member of each peri- 
odic orbit family encounters the singular vertex. Further, a non-periodic orbit originating 
from the same point but with a momentum slightly different from a periodic orbit suffers a 
net transverse deviation that equals (— 1)™^ sin(0 — P )Z</,. Here 1$ is the distance traversed 
by a non-periodic orbit at an angle after reflections from the boundary and (p p is the 
angle at which a periodic orbit exists. These facts can be used to converge on periodic orbits 
rapidly and the method works for other polygonal billiards. Details of this and the other 



method employed can be found in [17 



Note that in both cases, the curves in Fig. 2 are linear as expected from eq. (^j. For 
the integrable case, b = 0.25, while for the pseudointegrable billiard b ~ 0.27 |TB|| . This 
is the first difference between the two cases. The higher value for the pseudointegrable 
billiard is possibly due to the fact that at the singular vertex, there can exist more than one 
periodic orbit with the same value of <p p . It could also reflect diffraction effects that have 
been neglected in Eq. (||). These issues will be discussed in a future publication [jP^[ . 

It is clear then that the leading term for the counting function, N(l) is nbol 2 / (a(l)) 
with corrections provided by the the quantum energy eigenvalues. This result holds for all 
rational polygons including those which are neither integrable nor almost-integrable. Here 
the average projected phase space area, (a(l)) = (J2 P J2r a p)/N(l) where the summation 
extends over all orbits with rl p < I. For rectangular billiards, (a(l)) ~ AA and this gives the 
quadratic law for N(l) PJT6[]. For pseudointegrable billiards, we plot (a(l)) in Fig. 3. The 
saturation for large, I implies N(l) ~ I 2 asymptotically For smaller values of I, (a(l)) 
increases, indicating a (local) sub-quadratic law for N(l) that we have indeed verified. 

Note that the maximum area, a max , occupied by a family is four times the area of the 



billiard since the invariant surface consists of four sheets corresponding to the four possible 
momenta a trajectory can have. The value at which (a(l)) saturates is thus far smaller than 
a max . Hence the density of orbit lengths for the pseudointegrable billiard is far in excess of 
an equivalent integrable billiard having the same area. 

With these findings on periodic orbits, we now turn to the statistical analysis of the 



quantum spectra. A commonly used measure is the spectral rigidity, A(L) defined as [[19 



fj pL/2d av 

A(L) = (min ^ / [N{E + E) - a - bEfdE) (4) 
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where N(E) counts the number of eigenvalues, E is the energy at which the measure 
is evaluated and (.) is an averaging in energy over scales larger than the outer scale |T5 



determined by the slowest frequency of oscillation in Eq. ([I]). It is possible to analyze the 
rigidity in terms of periodic orbits using eq. ([I]) and the basic semiclassical expression for 
the rigidity is then JL5J : 



A(L) = <E E 4^ «s(fii - S 3 )H l3 ) (5) 

where the summations extend over all periodic orbits @, A = (a J 2 /327r 3 / iV ^) 1 / 2 ,T, = 
OSi/dE evaluated at Eq and Si = y/Eok. The function H i3 = F(yi — yj) — F(yi)F(yj) — 
3F'(yi)F'(yj) where yi = LTi/2d av and F(y) = sm(y)/y. 



For L 2nd av /T min and large E , Eq. (§) can be simplified further to yield [^5 
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A(L) = — / o -K{r)G{,Lr) (6) 

where G{y) = 1 - F\y) - 3F' 2 (y), r = T/(2nd av ), K{t) = 2tt0(T)/4, and 

<P(T) = (E E M, cos(^ - S 3 )5(T - (T t + Tj)/2)) (7) 
* i 

Eq. (H) is useful for analytical studies only when the collective properties of periodic orbits 
as embodied in 4>(T) are known. A first step in this direction is the diagonal approxima- 
tion [T5j for small T based on the fact that orbit pairs have large action differences and hence 



off-diagonal terms do not survive averaging. The diagonal sum, 4>d{T) = (J2iA?5(T — TA) 
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is however unknown for generic pseudointegrable systems, though on treating the area, dj, 
as constant and using the asymptotic law, N(l) ~ I 2 , one might infer that 4>d{T) is constant 



H as in integrable systems fl5|,|6|j. These assumptions are however incorrect especially for 
values of T where the diagonal approximation is expected to be valid. We shall therefore 
investigate this numerically 

For large T, off-diagonal contributions are generally important and more difficult to 



estimate ||21|| . For integrable systems where the density of orbit lengths is small, off-diagonal 



contributions vanish even for large T ||15|| . One might expect this to be true even in the 
pseudointegrable situation due to similarities in the asymptotic proliferation laws though 
it must be noted that the density in the pseudointegrable case is larger compared to an 
equivalent integrable system having the same area. 

In Fig. 4, we plot the function, J(r) = J T Kij'^dr' (upper curve) as well the diagonal 
part, Id(j) = Jq Kd(t')(1t'. The fact that Id(j) and I{r) coincide until r c ~ 0.42 implies 
that off-diagonal terms do not contribute for r < r c . The diagonal approximation thus 



provides a good estimate of K(t) for short times fZ2 



The important departure (as far as A(L) is concerned) from integrable behaviour lies in 
the fact that Id{ t ) displays variations in slope for r < r c indicating that K£>(t) is not a 
constant p3| . For r > r c , Kd(t) is constant but has a value much smaller than unity (for 



integrable systems Kd(t) = 1 for all values of r). Of significance as well is non- vanishing 
contribution of the off-diagonal part (/(t) — Id(t)) for r > r c . This is probably the first 
example where a power law proliferation of periodic orbits gives rise to non-zero off-diagonal 
contributions in the form factor. 

Finally, we compute the spectral rigidity using Eq. (]5|) with 6621 periodic orbits and 
estimate the contribution of the longer ones using Eq. (0) and an interpolation for K(t). 
The result is displayed in Fig. 5 where we also plot A(L) obtained numerically using the 
quantum eigenvalues. The agreement is remarkably good leading to the conclusion that 
diffraction effects are not as significant for the statistical properties of the spectrum. 

In summary, we have brought to light several interesting properties of periodic orbits 



in pseudointegrable billiards and numerically established that periodic orbits provide good 
estimates of spectral correlations even when diffraction plays a role. Our specific conclusions 
are the following : 

• orbit families obey the sum rule (X^>X^li a P S(l — rl p ) / \rl p )) = 2nb thereby giving rise 
to the proliferation law N(l) = 7r6 / 2 / (a(l)) for all rational polygons. 

• (a(l)) increases initially before saturating to a value much smaller than a max ; the 
asymptotic proliferation law is thus quadratic even for systems that are not almost-integrable 
and the density of periodic orbits lengths is far greater than an equivalent integrable systems 
having the same area. 

• the diagonal part of the form factor, K{r) approaches a constant at r = r c and the 
value is much smaller than unity. In contrast, the diagonal part in integrable systems is 
identically equal to 1. 

• off-diagonal contributions in K{r) are non-zero for r > r c even though the proliferation 
of periodic orbits is quadratic as in integrable billiards. 

• periodic orbits provide good estimates of correlations in the quantum spectrum. 

It is a pleasure to acknowledge useful discussions with Predrag Cvitanovic, Bertrand 
Georgeot, Gregor Tanner and Niall Whelan. 
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FIGURES 
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FIG. 1. A pseudointegrable billiard. Periodic orbit families are restricted in extent by the 
singular (37r/2) vertex. 
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FIG. 2. The function, S(l) = J2 P J2 r w f° r a rectangular and a pseudointegrable (upper curve) 



billiard. The sum is restricted to orbits with rl p < I. 



11 



<«(/)> 




0.7 1 1 1 1 1 1 

20 40 60 80 100 

I 



FIG. 3. The average area, {a(l)), as a function of orbit length. Note the saturation for long 
lengths. 
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FIG. 4. A plot of I(t) = Jq K(t')<1t' (upper curve) along with the diagonal part. Note the 
crossover at r c ~ 0.42 where off-diagonal terms begin to contribute. 
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FIG. 5. The spectral rigidity, A(L) for the L-shaped billiard. The lower curve is evaluated using 
periodic orbits while the diamonds are the values computed directly using energy eigenvalues. The 
straight line above is the Poisson result (L/15). 
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